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I.   INTRODUCTION 

In  this  piper  we  shall  study  strong  solutions  to  the  Cauchy 
problem  for  n  certain  class  of  non-linear  partial  differential 
equations.   In  genenl  the  problems  we  consider  pve   improper  in 
the  sense  of  Hadamard    .   That  is  the  solutions  do  not  depend 
continuously  on  the  initial  data  in  any  reasonably  defined  norm. 
The  general  class  of  problems  treated  includes  as  special  cases 
such  classical  instances  of  improperness  as  the  Cauchy  problem 
for  elliptic  equations  and  timelike  initial  value  problems  for 
hyperbolic  equations. 

There  are  at  least  two  ways  of  approaching  improper  problems 
T-dth  the  aim  of  obtaining  solutions  which  do  possess  satisfactory 
continuous  dependence  properties.   One  is  to  seek  solutions  only 
within  a  restricted  class  of  functions  (e.g.  functions  with  a 
fixed  finite  bound  in  a  suitably  chosen  norm).  This  approach  has 

been  extensively  investigated  by  P.  John  ^'nd  quite  complete  resulti 

f2l 

have  been  obtained    , 

The  second  method  which  has  been  employed  by  Miranker  ^     and 
the  present  author'-^''  is  b-^sed  on  the  restriction  of  the  class 
of  allowable  initial  and  boundary  data.   In  the  latter  work  the 
Cauchy  data  is  assumed  to  be  band-limited,  that  is  to  possess 
Fourier  transforms  with  compact  support  in  all  variables.   This 
assumption  is  motivated  by  consideration  of  the  limitations  of 
measuring  systems  which  provide  the  initial  or  boundary  data  in 
"lany  actual  physical  problems.   Under  such  a  stringent  assumption 
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it  is  not  very  surprising  that  very  strong  results  can  be  obtained. 
In  fact  we  have  shc-m  in  [I|-]  for  a  cl^ss  of  nonlinear  problems 
analogous  to  the  ones  treated  in  this  paper  that  the  resulting 
continuous  dependence  of  the  solution  is  of  the  Lipschitz  type. 

Several  questions  arise  naturally  in  connection  with  the 
method  of  restricting  the  class  of  allowgble  data:  (l)  when  the 
operator  involved  in  the  differential  equation  is  of  a  specified 
type  (e.g.  hyperbolic),  can  the  assumption  of  bandlimitedness  be 
limited  to  a  subset  of  the  independent  variables  and  strong 
continuous  dependence  of  the  solution  still  be  demonstrated;  (2) 
more  generally,  when  the  type  of  the  operator  is  unspecified,  can 
bandlimitedness  be  replaced  by  a  weaker  assumption  on  the  behavior 
of  the  initial  data  without  losing  the  Lipschitz  continuous  depend- 
ence of  the  solution. 

In  [5]  ve  have  given  a  partial  answer  to  the  first  question 
by  showing  that  for  timelike  initial  value  problems  for  hyperbolic 
equations  it  is  sufficient  to  assume  that  the  initial  data  is 
bandlimited  only  with  respect  to  the  space  variables.   The  result- 
ing continuous  dependence  of  the  solution  is  still  of  the  Lipschitz 
type. 

In  [6]  we  continued  our  study  of  timelike  hyperbolic  problems 
by  combining  the  results  of  [5]  with  a  relaxation  of  the  assumption 
of  bandlimitedness.   Instead  it  was  assumed  that  the  Fourier 
transforms  of  the  initial  data  "ith  respect  to  the  space  variables 
possessed  a  certain  exponential  growth  at  infinity.   Moreover,  we 
showed  in  [6]  that  this  specified  growth  was  the  best  requirement 
of  this  type  that  could  be  placed  on  the  initial  data  and  still 
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obtain  solutions  with  Lipschitz  continuous  dependence. 

In  this  paper  ve  shall  treat  a  class  of  nonlinear  (more 
precisely,  quasi-linear)  Cauchy  problems  which  are  similar  to  those 
discussed  in  [tj.],  and  in  T'hich  the  operator  type  is  unspecified. 
^'e  replace  the  bandlimited  assuraption  of  [[|.]  by  a  growth  condition 
on  the  Fourier  transforms  of  the  initial  data  (with  respect  to  all 
variables)  analogous  to  that  eiTiployed  in  [6].   Our  present  results 
together  with  those  of  [6]  provide  a  fairly  general  answer  to  the 
second  questione  raised  above. 

The  writer  is  indebted  to  Professors  F.  John  and  J.  Moser  for 
several  helpful  discussions  and  suggestions. 

II.   PORMULATIOR  OF  THE  PROBLEM 
Let  k  be  an  integer,  a  a  positive  real  number  and  <^   a  real 
variable.  Vb   define  a  function  ^,{   '—'  )  by 


'k 

a|uL,.,k 


(1)   4k(  ^^  ) 


,  W^   >  1 


CV   <  1 


Note  that  4^(  ^^  )  is  continuous  for  -co  <  ^J  <  oo  .  He  say  that 

a  real  valued  function  g  of  a  real  variable  x  belongs  to  the  class 

eJ(x)  if 

(a)  g  ^  L2(-oo  ,00  ) 

(b)  The  Fourier  transform  G(  CaJ  )  of  g(x)  satisfies 

G(  CJ  )  4^(C.;  )  €  L2(-co,co  )  . 

From  (b)  we  hnve 

(b')   There  exists  an  l^^   function  gj^(x)  such  that 

G(  U    )  4^^(  '^   )  =  G^(  (,J  ) 


•^^cf; 


where  G.  (  ^-^  )    is   the   Fourier  transform  of  g,  (x). 
Explicitly, 

r    e^l^^^l    G(ej)c.J^      ,     h^l    >1 

\     e^G(  CJ  )  ,    |tJ|    <  1      . 

From  these  conditions  it  follows  that  g,  when  thought  of  as  a 
function  of  the  complex  variable  z,  is  analytic  for  I Imag  z|  <  a. 
The  class  of  transforms  of  functions  of  E  (x)  will  be  denoted  by 
C-  ,  (  ^'  ) •   A  real  valued  function  g(x,y)  of  two  real  variables 
x,y,  where  y  is  restricted  to  the  interval  0  <  y  <  h  will  be  said 
to  be  of  class  E^  j^(x,h)  if  for  every  y  <=:  [0,h]  g(x,y)  e  E^(x), 

g(x,y)  z:  C^ij),    and  if  the  Lp  norm  of  g  with  respect  to  x  is 

CD  2 

uniformly  bounded  for  y  e  [0,h]  (i.e;,  ||g||  =  sup    (   lg(x,y)|  dx 

=  sup    I  |g(y)  I  I  <  CO  ) .    It  is  clepr  that  E,  „(x,h)  (and  its 
~  0<y<h  ^'"^ 

class  of  x-Fourier  transforms  C  \r  m   ^  '^  »h))  is  closed  in  L^  , 

for  each  fixed  y  e:  [0,h]« 

For  an  arbitrary  L^ ( 

transform  H(  CJ  ),  i-je  denote  by  [h(x)]-  (P  positive  and  real)  the 

inverse  transform  of  the  function  H(  C<-^  )e""  '   '.   Thus  [h(x)]o 

is  a  damped  or  smoothed  version  of  h(x).   This  operation  may  be 

thought  of  as  a  convolution  in  x-space 

h(x)  -:c-  d(x) 

where  the  smoothing  function  d(x)  is  given  by 

2(3 


For  an  arbitrary  L^(-co,co)  function  h(x)  with  Fourier 


d(x) 


«2  2 
p  +x 


Let  the  linear  operator  L  be  defined  by 


2l 


p+q 


(2)  L=i:    a   (y) 

P-qfm       ^^p^^q 

where  ^^Aj)  ^   C  (y)  for  0  <  y  <  h  and  a„   =  1.   Our  goal  is  to 
prove  the  existence,  uniqueness  and  continuous  dependence  of  L„ 
solutions  to  the  equation 

(3)  Lu  =  [f(u)]p 

in  an  infinite   strip   R(h    )   defined  by-oo    <x<od,    0<y<h, 
given  Cauchy  data 

(1;)  A!liL|z£l     =  V    (x)         ,      q  =   0,l,...,m-l 

on  the  X-axis  satisfying  v  (x)  e  ^m-a^"^^  "^°^  ^^"^   ^'   ^®^®  °-  ^^^ 
P  are  positive  real  constants        (to  be  determined  later) 

T'hich  depend  on  the  coefficients  a   ,  and  [f(u)]o  means  that  f(u) 

pq  P 

is  damped  ^s  a  function  of  x. 

It  lA'ill  be  assumed  th.?t  f(u)  satisfies  the  following  condi- 
tions (which  are  not  all  independent)  as  a  function  of  u: 

(1)   f  €  C°(u) 

(ii)   f(u)  has  the  quadratic  property  in  a  sphere  S^  of 
radius  \  (to  be  explicitly  determined  later).   That  is,  if 
I |u| I  <  \,  then 

|f(u)|  <  A^u-" 

where  the  constant  A,  depends  only  on  the  sphere  S   . 


This  condition  is  essentially  equivalent  to  assuming  that 
f(u)  <^  C^(u)  -  see  Ref  [!;]• 


(iii)   f(u)  satisfies  a  Lipschitz  condition  in  S,  ;  i.e;, 

A 

if  u,u'  e  S^,  then 

|f(u)  -  f(u') I  <  K^  |u-u'|  , 

where  K,  depends  only  on  S,  #    It  follows  immediately  from  (ii) 
that  f(0)  =  0  gnd 

00  00 

J   |f(u(x,y))|dx  <  A^  I   |u(x,y)|^dx  =  A^||u(y)||^  . 

-OD  -00 

Thus    if  u  <r  S- ,    the  Fourier  transform 
A. 

oo 

(5)  Pu^^'2r)=|^     j      f(u(x,y))e"^'^''  dx 

-oo 

exists  and  [f(u)]-  may  be  explicitly  defined  by 

-00 

We  denote  the  transform  of  [f(u)]-  by  [F  (  ^^  »^^^S  '  ^'^^^   (5) 
it  follows  that  P  (  C<J   ,y)  satisfies  the  inequality 

(7)        |f^(cj  ,y)|  <  2I-  A^||u(y)||2    . 

Ill   AN  EC^UIVALEWT  INTEGR/^L  EQUATION 

Our  proof  of  the  existence,  uniqueness  and  continuous  depend- 
ence of  solutions  to  (3),  (I4.)  will  depend  upon  the  formulation  of 
an  equivalent  integral  equation. 

First  we  formally  write  down  the  reduced  (ordinary)  differen- 
tial equation  dissociated  with  (3),  ik) 


*YKejiw 


where 


m-q 


(9)  o^.,(t^,y)  =Z:^  ap^(y)(i..-)P 

•1     CD  J 

(10)  U(uv/,y)  =  ^   (   u(x,y)e"  '^'  ""   dx 


-00 


and 


(11)   ^!u(_.oj_  =v  (lJ)  = 


"'  y 


2i 


v^(x)e"^  ""  ^  dx 


for  q  =  0,1,...,  m-1.   Note  that  c.  (  (,--'  ,y)  is  a  polynomial  in  i,J 
of  degree  k. 

Our  program  will  be  to  use  the  variation  of  parameters  method 
to  obtain  9n  explicit  integral  equation  for  U  the  solution  to  (8), 
(11).   To  do  this  we  first  consider  the  solution  Uq(  '-  ,y)  to 
the  homogeneous  equation 


(12) 


^  Uq  =  0   , 


SJ   u.,0) 


?.y^ 


V^(c^)   ;  q  =  0,l,...,m-l, 


Equations  (12)  can  be  written  as  a  system  and  estimates  of  the 
growth  of  VA  tJ>  ,y)  obtained  in  the  following  way.  ^/Fe  distinguish 
two  cases:  Case  A;  |   I  >  1.  Define  the  vectors  ^-^q  ( --^ ,  y ) ,  V(  C-'  ) 
and  the  matrix  ^'\i\   (  l-^  ,y)  by  _ 

Uo(  ^  ,1) 
1   3UQ(u^,y) 


(13)  ^'n^  '-^y)  = 


,ra-l 


-\  „m-l 


,  V(u;)  = 


Vo((J) 

i  V-(cJ) 

1-'  -■■ 


-i-T-  V  Au) 
m-1  m-1 


(li+)   !\(\   (  tJ  ,y) 


0       10 
0        0         1 


c  c      _ 

m         m-1 

,m-l 


^^ 


1 


-  '"l3^ 


Then   (12)    c?5n  be  written  in  the   form 

(15)  _z= =  tJ     •;|(w,y)  w  (eJ,y)    ,     ^■'„(CJ,0)  =  V(c^) 

dy  _2  _u 

or  equivalently  ^-^-Ci^^y)  can  be  expressed  as  the  product 

(16)  ^(  UJ  ,j)    =  Ur  (  L^,y)  n-.J)      , 

where  the  m?trlx  IV  (  CJ  ,y)  satisfies  the  rmtrix  differential 
equation 

(17)  ^^^-^^^-^     =    c^%   (  c^  ,y)  U^   (  (V  ,y) 

dy 

and  1^^  (  Cj   ,0)  =  ,J  ,  the  identity  matrix  of  order  m. 

Since  c,  is  a  polynomial  of  degree  k  in  i-^'  and.  is  bounded 
in  y  for  0  <  y  <  h,  we  have  that  It]  is  bounded  uniformly  for 
I  IV  I  >  1   and  0  <  y  <  h  •  Call  this  bound  |j,.  An  upper  bound 
for  ij,  in  terms  of  the  coefficients  a   can  be  obtained  explicitly: 


■:T0 


''\  i' 


o         m  o        o    m    c    o 

i,J=l   ^J  j=l   c-^ 


o   m  o ..  i 

=  (m-1)'^  +  E  (,.;  J  I  1_ 

0=1  _  p=o 

<  (m.l)2  +f^  [  f:  la      |] 

7^1  p=o  P'^'-J 


2 


,-(y)(i-)P|^ 


\''e   use  the  same  symbol  |  |  |  |  for  both  the  m-^trix  and  Lp  norms 
It  will  always  be  clear  from  the  context  which  is  meant. 


where  a  „  ^  =  max   a^  ^  ,-(y)»   Then  from  (17)  and  repeated 
P'"^-J   0<y<h  P'"'"^ 

application  of  the  Cauchy  inequality  it  follows  that 

(18)    II  ]v'(c..,y)n  <  II  3  iie^^yK-i  =^^e^7\'-\ 


so   that   from   (16) 

(l8a)      lh'^(^-.',y)ll  <   II  l-''(L^,y)IM|v(t^')M 

M7 


<y.e-'^-^'|r        c.-2Mv^(.^)l^]    ^ 


2_ 

q=o 


In  particular  then  for    |w  I    >  1 
(19)       |UQ(c-,y)|    <ym  e^'3^l'-'' 


q=o  "1 


or  xoniformly   for  0  <  y  <  h 

(20)         |u^(>^,y)l    </m  el^h'^^l    j   £'   c.-2<l|v   ( u' )  l^j 


m^l 
q=o 


Case  B;   j'-''  I  <  1;  here  T,.e  write  (12)  as  a  system  in  a  slightly 
different  way.   Define  the  vectors  U„  ( •'---^ ,  y ) ,  V(  ^-^  )  and  the  matrix 
'\\   (ej,y)  by 


(21)   UQ(u;,y)  = 


Uo(u>,y) 


Uq(  '^,y) 


ay 


.^-\(^,y) 


^y 


ra-l 


,  V(U7  ) 


m-i 


(22)      "Y^  (LJ,y)    - 


0        1        0 0 

0        0        1    0 

0        0        0 1 

c        c      . c. 

m  m-1                          1 
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Then   (12)    takes   the   form 

^^o^-^'^y) 

(23)     -^^ =     it   C^y)  u^(c^,y)     ,     u^(,'J,o)  =  v(L-) 

or  equivalently 

(2U)     Vq{^'^,j)   =    U  (^^,y)  V(c^) 

where  the  matrix  'V  ((^,y)  satisfies  the  differential  equation 
(25)        ^  ^^^  ^"^y)    =   n  (,.,y)  U  (W,y)  ,   U  (c^,0)  =  3  . 

dy 

Since  l'^-*  I  <  1,  the  matrix  M  is  uniformly  bounded  for  0  <  y  <  h. 
In  fact  it  is  easy  to  shoi'  that  n"\  and   \  have  the  same  bound. 
Thus  proceeding  as  before  ^'e  find 

m-1 

^■ 


(26)      l|UQ('.-^,y)M  <^/me^'y  \f^    |v  (^^)|^] 


q=o 


or   in  particular  uniformly  for   0  <  y  <  h 


m-1 


1 
2  ]     2 


(27)  \vJi^',7)\  l^^o^    \JZ     |v(tJ)rj 

^  *-  q=o        ^ 

Recalling  the  definition  (1)  for  (|)^(^)  we  can  combine  the 
inequalities  (20)  and  (2?)  by  defining 

(28)  a  =  tih 


11 


and  writing 


(29)    (UQ(^-^y)l  <  MUQ(t^,y)||  <ym]2_   l4!-( '^)V  ( ^^  )  T  , 
—  q=o    4     4.      J 


1 
2>  2 


The  solution  U  (0,y)  to  the  homogeneous  equation  (12)  is 
given  either  by  the  first  component  of  (16)  which  is  actually 
valid  for  all  y  ^'nd  any  iO  -region  excluding  a  neighborhood  of 
^•^  =  0,  or  by  (2I|.)  vhich  is  valid  for  all  y  and  «.^  •  The  purpose 
for  the  formal  separation  into  the  Cases  A  and  B  is  to  obtain  the 
estimate  (29)  which  is  sharper  than  would  result  from  a  considera- 
tion of  the  system  {2l\.),    (25)  alone.  A  similar  procedure  will  be 
employed  in  considering  the  solution  U((-^,y)  to  the  nonliomogeneous 
problem  (8),  (11). 

Define  the  vectors  U('^-^,y),  I:l(c^,y)  and  t  ^^^ '''^'^^  ^3  ^^ 
U(uv,y) 

1  ^  u(^,y) 

^    >7 


(30)   V(C^,y) 


1 

.  ,m-l 


u(^^'.y) 


:m=i 


U  (U,y)  = 


U(tJ,y) 
au(^,y) 

a  y 


?^"^"^U(CJ.y) 


3y 


m-l 


[P^O^lp 


[F,(u,y)], 


Then  (8)  may  be  written  as  a  system  in  the  following  way 


^i^ 


^■:;oJ:^;ol 


12 

(31)    =  ^-^    ^^^\   (^^,y)  v(^,y)+[F  (<.^,y)]  ,  H(cJ,o)=v(c-'); 

dy  -^ ^ 


du(L^»y)       v^ 

(32)      =    \|    (i^,y)   U('-J,y)+[P   (iV,y)].,   U(:y,0)   =  V((-.^'); 


|tJ|    <  1    . 
Consider  the    integral   equations 

y 

(33)      i:(<v,y)    =  WQ(c..,y)    +J      -lvMu,y)    Ll,''"^(i^»    )[F^(--%n  )3pdvi 

0 

.y 
(31!-)  u('<^,y)  =UQ(c-,y)   +   j     U(..vJ,y)  U"^(  .^  ■■V[F^(^---, -y  ]pd  r|   , 

where  U^  and  W  are  defined  as  before.   Prom  the  classical  theory 
of  ordinary  differential  equations  U_,  ^'q ,  "lA  and  lO^ ^  C      (y),  so 
we  may  different i'^te  in  (33)  snd  (3l4-)»   It  folloi^.'s  immediately 
that  these  expressions  satisfy  (31)  and  (32),  respectively,  includ- 
ing the  required  Initial  conditions.   As  above  we  note  that  the 
representation  (33)  is  valid  in  any  region  vhich  excludes  a 
neighborhood  of  Co  =  0,  and  (3i|)  holds  everywhere. 

If  we  denote  the  elements  of  the  matrices  'A  and  Ik        by 
U^  .  and  U  "1^  respectively,  then  using  {3k)   the  first  component 
'^i'-J  ,y)  of  U  may  be  expressed  in  the  form 

/^  m  -. 

(35)  U(Co,y)  =  UQ(u',y)  +     Z.     Uij.(L^,y)  uyc.,v-))[F^(..,^')))p 

Now  let 

00 

Uo(x,y)  =  \  Uo(C.,y)  e^'^'^  d  ^ 

-OD 


CD  13 

u(x,y)  =  I    U(LJ,y)e   ^  d  u^ 

CO 


u^j(x,y)  =  J   U^j.(i.J,7)e"^  ^  dc^ 


CO  ;■  j  =  1,2,  ...,m 


"li  <='•'''  =  J   V'^-y'^'"''''* 


•00 
CD 

[f(u(x,y))]p  =  j    [P^('-.',i^)]3e^''^  d.-.  . 

Then,  assuming  that  p  can  be  chosen  such  that  the  integral  in 

Ok)    (and  therefore  in  (35))  is  finite,  it  follows  from  (35)  and 

an  appeal  to  Pubini's  theorem  that  u(x,y)  satisfies  the  integral 

equation 

y  OD 

(36)  u(x,y)  =  UQ(x,y)  +(      1    k(x-^,y,  ■]  )[f  (u(4,  ^]  )  ]pd^d  v^  , 

0  -  00 

where  k(x,y,  "A  )    is  the  sum  of  convolutions: 

(37)  k(x,y,  n  )  =  -^ ^     u  (x,y)  ^- u";;  (x,  \)   )   . 

Thus,  when  p  is  chosen  in  this  manner,  u(x,y)  defined  by  (36) 
satisfies  the  original  equation  (3)  ^nd  the  boundary  conditions 
(1+),  since  U(  ^  ,y)  as  given  by  {3^)    satisfies  (8)  and  (11).   ^'e 
have,  therefore,  reduced  the  origin^^l  problem  to  proving  the 
existence,  uniqueness  and  continuous  dependence  of  a  solution  to 
the  pair  of  integral  equations  {3^),    (36),  or  more  gener?illy  to 
the  pair  (31;),  (36).   This  proof  is  given  in  the  following  section, 

IV   SOLUTIOI'I  OF  THE  INTEGRAL  EOUATIOWS 

^'fe  consider  the  integral  equations  {3k)  f    (36).   For  p  =  0,1,.. 

define  the  iterates  U  ,  u  by 
_P   P 


(38)  u  -^(^^,y)  =v^{u,j)  +   J  U  ((-,y)  U'^(v.,v^) 


Ik 


[F^(c^.',\^)]pdn 


y   00 


^^"^^      ^p+l^^'^^  =  UQ(x,y)  +  J   J    k(x-C,y,  v-)  )[f(Up(?,  n)]pd^dO 


0   -CD 


where 


Up,,(  ^-,y) 


Up^,(C^,y) 

?y 


m-l 


u„u.T('-^»y) 


a     "^p+i 


^y 


m-l 


([1.0)   Up^^(CJ,y)  =  UQ(CJ,y)  +  <|   ^  U^  j(CJ,y)U-^(  C^  ,  ^  )  • 

"o  J~^ 

In  order  to  obtain  sharper  estimates  in  whp.t  follOT-'s  ve    shall  use 
the  representation  (38)  for  |'-^|  <  1  and  for  \i-^'\   >  1  the 
equivalent  form 

(1^1)  ^_^(c^,y)  =  V'^'^^  ""  j  ^^^  ^^'"'^^  \^'"^(^-,  n)  . 

0 

_-E ^ 


The  expres3ion5(38)  and  (kl)  are  equivalent  for  |(v.'|  >  1 
(actually  in  any  region  including  a  neighborhood  of  c^  =  0) 
since  U  ,  '-'  ,  '^vA,  ,  "\^^   depend  only  on  the  solution  to  the 

homogeneouT'problem* 
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where 


v^(c^,y) 


1  aVii^'^^ 


■^'\..^^'7) 


Jli. 


m-1 


Prom  the  properties  of  Uq,  ^''  ,  tA-  and  l^^'   we  see  that  U 

and  ^'  ,  satisfy  the  initial  conditions  (11)  and  u  ,  satisfies 
p+1  p+1 

(i|). 

We  next   obtain  bounds  for    Mil    (^^,y)    l/l"^(<--^,  ^^  )  1 1 

II  V  (^^y)  IJ^-^C^^,  V))||  . 

Case  A;    Consider  the  adjoint  equation  to  (17) 


and 


dV-V^^.y,  V  =  .  l.r^-((.),y,  i^  )  u  Tn  ( w  ,  ;^  )   ,  0  <  n 


•.  ^-1, 


It  is  easy  to  show  that 

lt^'''(  ^-^  ,y."^ )  =  V  ( w  ,  v/] )  'l^;""  (  C^' ,  l^  ) 
Therefore,  using  the  results  of  Section  III  we  have  for  I'jj  >  1 
that 

II  W'  ((.■-,y)-Ur"^L-',^)ll  !>  e^^y'^^' 

Prom  this  inequality  and  the  mean  value  theorem  it  follows  that 
(1|2)  ll^(U,y)||  <  l|WQ(c^,y)ll  +  y/m  e^l''hl[P^  ( (V,^(^^))]  J I  , 

where  ^  is  chosen  appropriately  depending  on  U  ,    0  <  ^  <  y  <  h, 
and  as  before  a  =  \ih*      But 
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p  p  '^ 

so   that  combining    (I4.2)    and    (l8a)    we  have 

(43)    il-     l^^^y^ll   l^^e^'"^'    ^\llu-^'^    |V   (u;)|2j2  ^ 


^i^  ■    q=0 

+y   l[F„  (^,^)]el] 
p  '^ 

Case  B;   For  (iv^l  <  l  ,  rie  proceed  similarly  to  obtain  the  estimate 

II  U  (CJ,y)U  '^u;,^)  II  </m  e^^y 

Thus  employing  the  mean  value  theorem  pind  the  results  of  Section  HI 
T?e  obtain  from  (38)  the  inequality 

(lU^)  IJU  .,(L>^,y)ll  <  >  e^  [^  J:  |v  (U)r  ^  +  y|[P,  (CJ,?)]J] 
-^ —  ^  q=o    4      ^  "p       ' 

"e  may  combine  (I1.3)  and  (I4-I4.)  by  using  (29).   One  obtains  for  the 

first  component  U  ,-,(^-^,y)  defined  by  (I4.O) 

ik^)     |Up^^(U,,y)|  <ym  [\^      U;q(C^)Vq(.o)|2j2  ^y3a-^(.^)^ 


q=o 


P 


where 

-TLdJ)  =  max  {  1,  hv|  j    . 
^•'e  have  assumed  that  v  ix)  ^   ^m-a^''^^  vith  a  given  by  (28). 
Thus  (|)^  V  ^  Lp  for  q  =  0,1,..., m-1  so  certainly  i^  V  e  L  for 

all  such  q.   Let 

(i;6)        C  =r   114  VJ  I 
q=o     "^  ^ 


Here  and  in  what  follows  the  double  b?r  denotes  the  Lp  norm  • 
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^^e  now  choose  P  to  be  a  real  constant  satisfying  3  >  a;  say 
P  =  a  +  5,  where  6  is  positive  but  otherwise  arbitrary.   Then 
from  the  definition  of  the  operator  [   ]_,  the  inequality  (7)  and 
the  Plancherel  theorem  we  have 

|[p^  (u;,^)]pl  <A^  IIUp(y)||V^I^I  . 

Therefore 


u^    P-     .'  .    -         '-Up 


ikl)  <  A^  IIU  (y.Mf  ;  -  , 


|L~^|<1  P 

I   I>1 

Combining  ikS) ,    (i|-6)  and  {l\-7)   we  obtain 

1 

MUp^l(y)||  </Sc  ^y/^A^    (  I  e-26  +  e^^)^  ||Up(y)||^ 

or  uniformly  for  y  e  [0,h] 

1 
(^-8)    MUp^JI    <ym  C   +  h/^     A^(   f  e-25   +  e^^)^    ||Up||2      ^ 


where 


IIU  II  =  max   ||U  (y)ll  • 
P     0<y<h   P 

Now  if  I |U  II  <  2C/m   ,  then  from  (U8)  it  follows  that 


|U  ^^M  <  2C/m  as  long  as 


^:^s. 


1  '' 

(!+9)  A^hP  2  (|  0-26  ^  g2aj  ij.c  <  l  . 

Thus  we  have  inductively  that  | |u  II  <  2C/ra  for  all  p  under 

the  condition  (l4.9)«   To  see  more  clearly  what  this  condition 

implies  for  h,  we  write  5  =  hv  (recall  a  =  hp,)  so  that  0  =  h(p,+v). 

Thus  if  we  explicitly  define  the  sphere  S.  ,  by  setting  X  =  /Sroir  C 

we  have  from  the  above  results  and  the  Plan'.cherel  theorem  that 

the  iterates  u  defined  by  (39)  all  lie  in  S,  as  long  as  h  is 
P  A. 

sufficiently  small  for 

(50)  v^h   (^  e-2h^  +  e^^^^)  <  ikCA^r^   /[T+v" 
to  hold". 

i''e  now  estimate  the  difference  ||u  ^,  -  U  ||  and  show  that 

p+1   p 

the  sequence  ■(  U  1  converges  in  Lp*   Proceeding  as  before  it 
follows  that 

p       '^     p-1     ^ 

(51)  =y/^e^'''^'-P''^'|F^(W,,)  -P^   (...5)1   . 

p       p-1 

Now  ail   -B|eJ|    =-5|lJ|    if    |uy|    >1   and   a  -  3  |  u/ 1    <  a  if 
1^  I    <  1.      In  any  event   e^       "^ '       '    <  e^  so   that   from   (5l) 

"Vl^^^    "  ^p^y^ll   ^  ^^^"^    ll^u   ^^^    -  ^u        ^^^"        • 

P  P-1 

Thus  applying  the  Lipschitz  condition  on  f  and  the  Plgncherel 
theorem  we  obtain 


Note  that  as  v  ->  oo  (in  which  case  (3)  degenerates  into  the 
linear  problem  Lu  =  0)  the  solution  given  by  our  method  reduces, 
as  it  should, to  the  unique  global  one  guaranteed  by  the  classical 
theory* 


ir-'^^rrr^;;/.:/ 


19 


"^p+l^y^  -  Up(y)||  <  h/m/27t  e""  ||f(Up(i:))  -  f  (Up^^(^^) )  |  | 

<  hv/i^72^  eV  llu  (^)  -  u   (^)ll 
-  A     p        p-1 

=  h/S  e''  K.  ||U  (S)  -  U  ,(^)|| 
\  p       p-i 

<h/^e<'K^  l|Up-Up.^|| 

for  all  u  e  S  ,  \  =  /Siun  C,  and, consequently. 

Therefore  if  (50)  holds  and  in  addition 

(52)      h  e^^"  <  [A  K^]"-^  , 

the  sequence  -\  U  \  defined  by  (38)  is  a  Cauchy  sequence  in  Lp  so 

that  the  limit  function 

U((J,y)  =  lim  U  (a;,y) 

p->CD  ^ 

(for  each  y) 
exists  belongs  to  L»[and  satisfies  the  integral  equation  (35) • 

Similarly  from  the  equivalence  of  (39)  and  (i;0),  and  the  uniqueness 

of  the  Fourier  transform,  the  limit  of  the  iterates  defined  by 

(39) 

u(x,y)  =  lim  u  (x,y) 
p->oo  ^ 

exists,  belongs  to  S,  ("L  )  and  satisfies  (36).   This  completes 

the  proof  of  existence  of  the  solutions  to  the  integral  equation 

(and  thus  to  the  original  problem  (3),  ik))    in  the  region  R(h), 

where  h  satisfies  (50)  and  (52).  There  remains  to  show  uniqueness, 

differentiability  (C  (x,y))  and  continuous  dependence  of  the 

solution. 
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The  uniqueness  of  the  solution  UQ(t^,y)  to  the  homogeneous 
equation  (12)  follows  from  the  classical  theory  of  ordinary 
differential  equations.   Thus,  to  prove  the  uniqueness  of  the 
solution  pair  lliCJ,j),   u(x,y)  it  is  sufficient  to  consider  a 
second  solution  U'((,J,y),  u'(x,y)  to  the  integral  equations 


U«(c^,y)  =Up{t^,y)  +  ^  ll(e>-',y)  XXj^^il^,^) 
0 

(53) 


[P^,(c^,v3)]pdr, 


y  OD 


u'(x,y)  =  UQ(x,y)  +J  \        k(x-§,y,  ^   )[f(u«(4,  -j  ))]pdCdP|   . 
0  -co 

Proceeding  as   above  we   find   from   (38)    and    (53) 

I  |U'  (y)    -  Up(y)  I  I    <  hj/iiT    K^e'^    |  |u'  (y)    -  Up^^(y)  I  I 

from  which  it  follows 

I  |U'  -  U  M  ->  0 

whenever  (50)  and  ($2)    are  satisfied.   Thus  from  the  uniqueness 
of  the  Fourier  transform 

||u'  -  Up  I  I  ->  0 

so  that  the  solution  u(x,y)  is  unique  in  Lp(-oo,+co),  yi[0,h],h 
satisfying  (50)  and  (52). 

It  follows  from  the  classical  theory  that  uA^^-^fj)    and 
lV  (UJ,y)  ^  C  (y).  Thus,  since  these  quantities  involve  derivat- 
ives of  U„(t-^,y)  of  order  m-1  in  y,  we  have  from  (35)  (or  (33)) 
that  U(LJ,y)  c  c'"(y)  and,  consequently,  u(x,y)  S   C  (y)  for 
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0  <  y  <  h.  To  prove  vl{x,j)   e  C  (x,y)  it  is  therefore  sufficient 
to  show  that 

i^k)  C)""  U(CJ,y)eL2   . 

Prom  the  definition  (1)  and  {l\.$)    ve  obtain 

ra-1  i 

Kj"^  U(u.--  ,y)  I  <  ym  j  C   'Cq^'-'  ^Vq^'-^'  ^  'l^  ^ 


(55) 


^-y/m  e^"^^^"^  '  t^u^'^' '^^  ^6  "  ^-^ 


Since  v  (x)  v^  ^m-n^"'^^  ^^  assiomption,  the  first  term  in  iSS) 
belongs  to  Lp  &nd  (5U)  will  therefore  follow  if  we  can  show 

(56)      ^aS^-{^)   ^,,m  |[p^(cJ,^)]J  eL^ 

for  all  Z   e  [0,h].   But 


and  thus 


lu'-'e"^  tV«"pll'i*X  ll°ll^   (■  /^''"e-2=l'-^U, 


|U|>1 


^        )        c^^e^to-Nt^l)  a^] 


uniformly  for  ^^  [0,h].   Since  the  integrals  are  finite  for 
5  >  0,  the  desired  result  (56)  follows. 

i''e  next  consider  the  continuous  dependence  of  the  solution. 
Let  U(/j  ,y),  u(x,y)  and  U(f'',y),  u(x,y)  be  solutions  to  (8)  and 
(3),  corresponding  to  the  initial  data  f  V  (u;  )  j  /^  ^  (x)  J   and 
iV  (tJ)  j  ,  j  V  (x)  ]  ,  respectively,   ve  assume  v  (x),  v  (x)  r 
e"^  (x)  for  q  =  0,1,..., m-1,  where  a  (which  depends  only  on  the 
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operator  L)  is  given  by  (28).   Thus 

q  =  0,1, . . . ,m-l.   Suppose  the  data  v  and  v  are  close  in  the 
sense  that 

(57)     nCq'V%)ii<e 

for   q  =   0,1, . . .,m-l,    and    some    c  >  0.      Then  certainly 

(58)  IU2q(Vq   -  Vq)||    <    e 

for  q  =  0,1,..., m-1.   Proceeding  in  the 

usual  manner  we  derive  the  inequality 

|u('J,y)  -U(cJ,y)|  <^m{i:   U""  J  C.)  )  1^  |v^(  OJ  )  .Y    {'J)\^j     ^ 

i     q_Q    -4  4  4        -^ 

+  y  /m  e^  "  ^"^'^  I  [P^(u;,^)]p-  [P-  (^-',^)]pl 
Thus  defining  the  sphere  S,  by 

\  =  max  [  1/8^  ZI  I  lOa"'  ^^  H    '  1 4^  J  J  '  ^ 


•q  q       ' —  ■  -q 


q"i 


and  assuming  that  condition  (iii)  on  f(u)  holds  for  u  t  S^,  we 
have 

l|U(y)  -  U(y)||  <\/m  j  £   •  '^""a^^a  "  ^a^  "^]  ^  +  y/m  e^  K^ 

'■   q=o     q  H    4 

||U(?)  -  u(^)ll  . 
^rhence  from  (58) 
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I  |U  -  ul  I  <  me  +  h/m  eV  I  |U-u|  | 


so  that,  if  (50)  and  (52)  hold, 

llu-ull 


1  -  h/m  e°7:^ 


Thus 


m/git _   m/2-R 


1-h/m  e^^       1-h/m  e^^^K^ 


and  the  solution  u  to  (3)  depends  Lipschitz  continuously  on  the 
initial  data  (1^.). 

Finally  vie   shov?  that  the  solutions  U  and  u  depend  continuously 
on  the  function  f(u).   Let  U(uJ,y)  u(x,y)  be  the  solutions  of  the 
equations  -4^  U  =  [G  (CJ,y)]   and  Lu  =  [g(u)]g,  respectively,  with 
the  same  initial  data  (11)  and  {l\.) ,     t'ithout  loss  of  generality 
ve  may  suppose  th^t  f  and  g  satisfy  conditions  (i),  (ii),  (iii) 
in  the  sphere  S.  ,  X  =  \/8Tim  C,  defined  by  (Lt.6).   Prom  condition 

A. 

(ii)  ve  have  that  f(u(x,y)),  g(u(x,y))  e  L.  as  functions  of  x 
for  u  S'  S.  .   '"'e  suppose  that  f,g  are  close  in  the  sense  that 

00 

(59)   max   i-  f    |f(u(x,y))  -  g(u(x,y))|dx  <  e 
0<y<h  '^^   ^^-00 

Then  using  the  same  estimation  procedures  as  before  it  follows 

that                      1               1 
v/2^^(vH-n)"  2(U  ^-2hv  ^  ^2htij2 
llu-ull  < ^ " 


1-hj/m  K^e^'^ 
as  long  as  (50)  and  (52)  hold. 


2h 

Ve  may  summarize  the  results  of  the  above  sections  in  the 
following 

Theorem;   Consider  the  Cauchy  problem  (3),  (i;),  vhere  the  operator 
L  Is  defined  by  (2).   Let  C  be  given  by  (I4.6)  and  the  sphere  S, 
defined  by  X  =  /8itm  C.   Suppose  |i,  is  a  bound  for  the  matrices  ''1 
and  l\  given  by  (II4.)  and  (22),  respectively. 

Then,  if  f(u)  satisfies  conditions  (i),  (ii),  (lii)  for 
ue  S^,  ^q^^)  ^  vcS^a      ^^  ~   0*l>««»»"i-l   ®^^  ^  satisfying  (50), 
(^2))  and  (3  is  chosen  larger  than  hjj,,  there  exists  a  solution 
u(x,y)  to   (3),(l|)  in  R(h)  =  -[  (x,y)  :  -00  <x<oo;  0<y_<hj 
vith  the  following  properties 

a)  u(x,y)  C  E°(x,h)  in  R(h) 

b)  u(x,y)  is  unique  in  the  class  of  functions  of  x  and  y 
which  are  in  Lp  with  respect  to  x  uniformly  in  y. 

c)  u(x,y)  depends  Lipschitz  continuously  (in  the  norm: 
/CO  2 

max  I  \         |u(x,y)|  dx]  ^)  on  the  initial  data  v„(x)  and  the 
0<y<h  --00  ^ 

function  f(u),  when  the  norms  of  these  quantities  are  defined  by 
(57)  and  (59),  respectively. 


The  extension  of  our  methods  to  the  Cauchy  problems  for 


the  equations  Lu  =  [f(u,^-^)]  ,  Lu  =  [f(u,  3^)]g  and  most  general- 


where  L  is  defined  by  (2),  can  be  carried  out  by  an  adaptation 
of  the  methods  of  Ref  [I;]  in  a  manner  completely  analogous  to 
our  present  treatment  of  the  problem  Lu  =  [f(u)]g  • 

J.M.  Zimmerman 
July  2,  1962 
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